Measurement of sub-shot-noise spatial correlations without subtraction of background 
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In this paper we present the measurement of sub-shot-noise spatial correlations without any 
subtraction of background, a result paving the way to realize sub-shot-noise imaging of weak objects. 
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Quantum properties of the optical field represent a 
resource of the utmost relevance for the development 
of quantum technologies, allowing unprecedented results 
in disciplines ranging from quantum information and 
metrology [T] to quantum imaging (for recent reviews see 
[2]). For applications in this latter field, a fundamental 
tool to be realized is the spatial correlation, at the level of 
quantum fiuctuations, of two optical beams. An example 
is the high-sensitivity detection of weak objects proposed 
in [2, which exploits the intrinsic multi-mode quantum 
correlations of parametric down conversion (PDC) emis- 
sion, a result that may have disruptive practical applica- 
tions in biomedical imaging or whenever there is the need 
of illuminating an object with a low flux of photons |20j . 
This technique is based on the spatial correlation in the 
quantum noise of two conjugated branches of PDC emis- 
sion IlllSlI^ . The image of a weak absorbing object in one 
branch, eventually previously hidden in the noise, can be 
restored by subtracting the strongly correlated spatial 
noise pattern measured in the other branch. In order 
to achieve a sensitivity superior than that available with 
classical techniques, one should be able to reach a sub- 
shot-noise (SSN) regime in spatial correlations even in 
the presence of the unavoidable background noise (elec- 
tronic noise, scattered light) [3]. Nevertheless, up to now, 
such a result was not yet achieved. Indeed, previous 
demonstrations of the quantum nature [H |9] of the spa- 
tial correlation of PDC beams were realized only for low 
photon numbers. In that regime the background noise 
was dominant, so that a proof-of-principle demonstra- 
tion of SSN correlation was possible only after correcting 
the results for the background noise (i.e. by subtracting 
the variance of the spatial pattern of the background, 
measured in the absence of PDC). Clearly, such a regime 
cannot be used for concrete imaging schemes, as the im- 
age distribution would remain hidden in the background 
noise. Similarly, a single-mode sub-shot-noise intensity 
correlation [lOl HH [HI US] [Ml [15] cannot be used to 
retrieve high-sensitivity information on the spatial dis- 
tribution of an object, since the quantum character of 
the correlation vanishes when one detects small portions 
of the beams instead of the whole beams. Thus, high 
sensitivity imaging requires spatial quantum correlations 
(i.e. SSN intensity correlations between several portions 
of two twin beams) in a regime where the photon flux 
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FIG. 1: Experimental setup. A triplicated Nd-Yag laser 
beam, after spatial filtering, produces type II PDC in a BBO 
crystal, which is then measured (after eliminating UV beam) 
by an high sensitivity CCD camera. 



is high enough to make the background noise negligible. 
Purpose of the present letter is to present this achieve- 
ment: a clear observation of SSN spatial correlation, 
without any correction for the background noise, is pre- 
sented and discussed, opening the concrete possibility of 
realizing imaging with a sensitivity beyond the standard 
quantum limit. 

The strategy we shall follow to reach our goal is that 
indicated by [3] : the major limitation of the experiment 
was the presence in each beam of a large excess noise 
with respect to coherent light, which made SSN corre- 
lation fragile with respect to unavoidable experimental 
imperfections (unbalances, errors in the determination of 
the symmetry center of the signal-idler distributions) and 
rapidly deteriorated the correlation for increasing gain. 
Such an excess noise can be lowered by working in a dif- 
ferent regime, i.e. a pump pulse of duration much longer 
than the typical coherence time of PDC beams, which is 
on the order of ps for a few mm crystal. As predicted 
in [3J, a nanosecond pump pulse should enable the ob- 
servation of SSN correlations at a much higher photon 
number. In our setup (Fig{l]) a type II BBO non-linear 
crystal {1 = 7 mm) is pumped by the third harmonic (355 
nm) of a Q-switched Nd:Yag laser (with pulses of 5 ns. 
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10 Hz repetition rate and a maximum energy of about 
200 mJ) after spatial filtering (a lens with a focal length 
of 50 cm and a diamond pin-hole, 250 /zm of diameter), 
for eliminating the non-gaussian components, and coUi- 
mation (by a system of lenses to a diameter of Wp — 1.25 
mm). After eliminating the pump (by a couple of UV 
mirrors, transparent to the visible, T=98% at 710 nm, 
and by a low frequency-pass filter, T=95% at 710 nm) 
the far-field of the down converted beams (signal and 
idler) is imaged onto a CCD camera, by means of a lens 
(/ = 10 cm) in a / — / optical configuration, ensuring 
that we image the Fourier transform of the crystal exit 
face. Therefore each transverse mode with wavevector q 
of the emission is associated to a point x = (A//27r)q of 
the image. We use a 1340X400 CCD array, Princeton 
Pixis400BR (pixel size of 20 /im), with high quantum 
efficiency (80%) and low noise in the read out process. 
The CCD exposure time is set to 90 ms, thus each im- 
age catched in this time window corresponds to the PDC 
emission generated by a single shot of the laser. After 
having identified the degenerate emissions (two separated 
rings in our phase matching configuration) , at wavelength 
\s = Xi = 710nm, by inserting an interference filter, we 
focus the analysis on two correlated areas of side 2,72 
mm within a 10 nm bandwidth around it. In each of 
the two signal-idler regions, the unwanted contribution 
of frequency components of the orthogonal polarizations 
other than the degenerate one are eliminated by a polar- 
izer T=97%). Notice that differently with the setup of 
[8] the signal and idler beams follow basically the same 
optical path, so that unbalances are reduced. 

In the ideal case of a plane wave pump (qp„„ip = 0), 
the momentum conservation in the microscopic PDC pro- 
cess requires that each pair of signal and idler photons is 
emitted in symmetric directions qs = — q^. As a conse- 
quence, perfect intensity spatial correlation between any 
symmetric portions — — x,j of the far-field zone of the 
signal and idler beams were predicted [H [51 [B]. For a 
finite pump of waist Wp, having a distribution of trans- 
verse momenta, there is an uncertainty in the relative 
directions of propagation of twin photons represented by 
the angular bandwidth of the pump, so that the photon 
numbers collected from symmetrical portions of the far- 
field zone are perfectly quantum correlated only when 
the detection areas are broader than a coherence area 
(roughly the transverse size of the mode in the far field) , 
whose size Aa;^^^ depends on the pump angular band- 
width \/wp [51 E], and on the parametric gain [511^ [IB]. 
Thus, we expect a quantum correlation in the number 
of photons Ni and Ng detected at any couple of sym- 
metric positions with respect to the center of symmetry, 
providing our detection areas are larger than the typical 
coherence area. 

The quantity that characterizes the amount of correla- 
tion is the noise reduction factor a, defined as the vari- 
ance of the photon number difference iV.^ — Ns normalized 



to the shot noise. In the ideal case of detection areas 
much larger than the coherence area we have the usual 
"twin beams" result: [11 [17] 



(1) 



where rj represents the total transmittance of the optical 
channel, including the quantum efficiency of the detec- 
tors. For classical light a is above 1 and reaches the low- 
est limit in the case of coherent light (cTcoh = 1) which 
is referred to as Shot Noise Level (SNL). Whereas for 
twin beams a is always smaller than unity (see EqjlJ. 
The other relevant quantity is the Fano factor Fg^i) = 
(^^iV5(i))/(7Vs(i)), which represents the fluctuation of the 
single beam normalized to the shot noise. For PDC , since 
each beam taken alone has a thermal-like statistics, the 
Fano factor has the multi-thermal form: 



l + E„ 



where = {Ng(^i))/M is called "excess noise" and M is 
the degeneracy factor, proportional to the number of spa- 
tial and temporal modes collected in the measurement. 

We estimate the expectation values as averages over 
a large spatial ensemble (large regions of a single im- 
age). Specifically, let us define iV(x) the number of pho- 
tons, registered by the CCD pixel in the position x of a 
generic region R of the CCD. SN{x) = N{x) - (iV(x)) 
is the fluctuation around the mean value, estimated as 
{N{x.)) = (l/jT-) X]x^(^)' with n the number of pixels 
in R. We estimate the a parameter by fixing a region Rg 
in the signal portion (136X136 pixels), and searching the 
corresponding region Ri of the idler field that minimizes 
the experimental noise reduction factor: 



,(C) = 



{S' [jV.(x)~iV.(-x + 0]) 
(iV,(x)+iV,(-x + 0) 



(2) 



where ^ is the displacement vector of Ri in the pixel space 
with respect to the optimal position. In other words, the 
numerator represents the residual spatial noise of the im- 
age in Rg after having subtracted pixel-by-pixel the cor- 
related noise pattern in Ri. aexpiO is plotted in Figj2] 
Far from the minimum, we are in the situation in which 
the two regions are completely uncorrelated, and aexp 
approaches the Fano factor of the single region (we sup- 
pose for simplicity Fg « Fi). Furthermore, Figj2] (a) 
shows that the size of the coherence area is larger than 
the pixel and the Full Width Half Maximum of the deep 
gives us Axcoh in terms of pixels. For this image, we es- 
timate Axcoh ~ 8 pixels (160/im). Since the SSN can be 
achieved only for detection areas containing one or more 
coherence areas, we performed a binning of the pixels, i.e. 
we group the physical pixels in blocks of dimension 8X8, 
called superpixels. The number of photons collected in 
a superpixel is the sum of the photons collected by each 
original pixel. Our CCD allows to perform the binning 
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FIG. 2: Example of noise reduction factor a^^p, evaluated 
for a specific image, as a function of the 2D position of the 
region Ri (left panels). The right panels are sections of the 
minimum, (a) a is measured using the CCD pixels (side 20 
fim). Its minimum value stays just 0.15 below the Fano factor 
{a evaluated far from the minimum) and is above 1. (b) a is 
measured using larger pixels (side 160 /xm). In this case, the 
correlation is higher (the gap between the Fano factor and the 
minimum is about 0.6) and, especially, the minimum is below 
1, i.e. the correlation is sub-shot-noise 



operation at the hardware level, also lowering strongly 
the signal to electronic noise ratio. Fig|2] (b) shows aexp 
in the case of binning. The deep aexjAQ), representing 
the estimate of a, is noticeably under the SNL while the 
FWHM of the coherence area is clearly on the order of 
the superpixel dimension. We note that, although the 
operation of binning reduces the spatial statistical en- 
semble, nevertheless the chosen region Rg is large enough 
to contain about 300 superpixels. This number provides 
an estimate of the spatial resolution, because it gives the 
number of details of an image that one would be able to 
distinguish in an imaging experiment. It is clear [5] that 
the resolution in this high-sensitivity scheme is limited 
by the size of the coherence area of PDC, which in prin- 
ciple can be reduced by acting on the pump transverse 
size (a systematic experimental study of this dependence 
having been presented in |16j). 

FigjS] presents the minimum values of Cexp plotted as 
a function of the photon number, where for each data 
point the center of symmetry of the two images (the po- 
sition of the minimum of aexp) is determined with sub- 
pixel resolution, using the algorithm described in [IS] . 
The corresponding Fano factor F — {Fg + Fi)/2 is also 
reported. Each point of the main graph represents the 
result evaluated for a single image, generated by single 
laser shot, while in the small insert the data are averaged 
over several images. The comparison between F and a 
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FIG. 3: The noise reduction factor aexp, without any back- 
ground correction, is plotted versus the mean photon number 
per pixel. Each point (crosses) in the main window corre- 
sponds to the analysis of a single image. The corresponding 
Fano factor, F, is reported by the circles. Several images reach 
a noise reduction under 0.7. The insert on the top-left angle 
presents the same data averaged over tens of images. The 
equation for the fit is y = (5.9 ± 1.5)10"^^ a; -f (0.78 ± 0.02). 



gives a measure of the reduction of the noise level that 
can be observed in a differential measurements, where the 
intensities pattern of correlated regions are properly sub- 
tracted. Although no a posteriori correction for experi- 
mental noise and imperfections has been applied, all the 
fitting line is clearly SSN in the range of photon number 
considered, with a minimum of 0.82 for a photon number 
per pixel of 600 (the smallest value for a single shot being 
0.66). Also the dispersion of the data is reduced, being 
almost all the point under the SNL. This figures are fur- 
ther reduced to cr = 0.7 when a 12x12 binning is used (the 
smallest value for a single shot being 0.5). This two facts 
represent a fundamental new feature with respect to [B] , 
where a lot of images were above SNL and a heavy noise 
correction had been applied on the experimental values 
of a. Our data demonstrate that in our setup the spatial 
fluctuations of two selected regions, comprehensive of the 
unavoidable electronic noise and diffuse light, are really 
correlated at the quantum level, not only in principle, 
and therefore it is possible to exploit this correlation for 
real and new application to imaging. For example for an 
object transmittance a < 15%, following the theoretical 
prediction of ^ we can estimate that our setup would 
provide a SNR more than 10% better than classical light 
with 8x8 binning (20 % with a 12x12 binning). 

Let us consider the sources of noise that can spoil the 
spatial correlation. First of all the electronic noise in 
the read-out and digitalization process is A = 4 pho- 
toelectron/pixel and A = 9 photoelectron/superpixel in 
the case of binning. Then, the quantum efficiency of 
the pixels fluctuates spatially with a standard deviation 
dt] = 3%. Anyway, this is a very minor effect in our pho- 
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FIG. 4: Noise reduction factor acor and Fano factor F, cor- 
rected for the background noise. Here we present a subset of 
the same data showed in Fig. |4] once the electronic and the 
diffuse hght noise has been subtracted. The equation for the 
fit is y = (3.7 ± 1.1)10"'' a; + (0.38 ± 0.08) 



ton number regime. On the other side, the fluorescence 
of the crystal, mirrors and filters, absorbing the high en- 
ergy pump pulse, together with the residual light of the 
pump itself, contribute to a background of diffuse light 
collected by the CCD in each shot. This represents the 
main source of noise in our present configuration. For 
the sake of comparison with previous experiments, we 
can correct for this background noise, acquiring an image 
where the PDC emission has been suppressed, simply ro- 
tating the pump polarization by 90 degrees. Let us define 
S'^j.^j the variance of the spatial fluctuations of this image 
in the regions Rs{i), and Mg^i^ the mean photon numbers. 



We have typical values of 5^^^-) 



200 and AL 



120. 



The formula for the corrected noise reduction factor is 



(iV,(x)-f iV,(-x))-M,-A/, 



(3) 



The same correction can be applied to the Fano factors. 
Figji] reports the results applied to a subset (for graph- 
ical reasons) of the data of Figjs] Some image reaches 
values of noise reduction of 0.4 and the average is un- 
der 0.6 in the lower photon number condition. The Fano 
factors are reduced as well. These data can be com- 
pared with the theoretical threshold for the noise reduc- 
tion given by Eq.([T]). Since we estimate in our setup a 
total transmission-detection efficiency of 67%, it leads to 
a = 0.33. This means that we still have a net margin of 
improvement of 0.27. We attribute the gap between the 
theoretical prediction and the experimental achievement 
to i) the fact that in our setup the superpixel is just on 
the order of the coherence area, while the optimal condi- 
tion for the observation of SSN would require a detection 
area sensibly larger than the coherence area, and ii) a 



residual error in the determination of the symmetry cen- 
ter of the signal-idler distributions, because the bigger is 
the superpixel, the larger is the uncertainty in the evalu- 
ation of the symmetry center, even when this is retrieved 
by using a sub-pixel algorithm, as described in [IB] . 

In conclusion, in this paper we presented our results 
showing that we reached the condition of sub-shot-noise 
in spatial correlation without any background subtrac- 
tion, a fundamental achievement for realizing high sensi- 
tivity quantum imaging protocols. Indeed, we obtained a 
substantial reduction of the noise under SNL in the spa- 
tial distribution of a large portion of PDC emission by 
subtracting the correlated one. 
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